The Rise and Fall of a Networked Society 
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We propose a simple model of the evolution of a social network which involves local search 
and volatility (random decay of links). The model captures the crucial role the network plays for 
information diffusion. This is responsible for a feedback loop which results in a first-order phase 
transition between a very sparse network regime and a highly-connected phase. Phase coexistence 
and hysteresis take place for intermediate value of parameters. We derive a mean-field theory 
which correctly reproduces this behavior, including the distribution of degree connectivity and the 
non-trivial clustering properties. 



Recent phenomenological studies on complex networks 
in social sciences have uncovered ubiquitous nontrivial 
statistical properties, such as scale free distribution of 
connectivity or small world phenomena Pj S Q ■ These 
properties have striking consequences on the processes 
which take place on such networks, such as percolation 
01, diffusion |^ |^, phase transitions and epidemic 
spreading 0- The research on complex networks raises 
questions of a new type as it addresses phenomena where 
the topology of interactions is part of the dynamic pro- 
cess. This contrasts with traditional statistical mechan- 
ics, where the topology of the interaction is fixed a priori 
by the topology of the embedding space. 

Phenomena of this type are quite common in social 
sciences where agents purposefully establish cooperative 
links Q. Links between individuals in a social network 
support not only the socioeconomic interactions that 
determine their payoffs, but it also carries information 
about the state of the network. This aspect has impor- 
tant consequences in the long run if the underlying envi- 
ronment is volatile. In this case, former choices tend to 
become obsolete and individuals must swiftly search for 
new opportunities to offset negative events. The role of 
the network for information diffusion is particularly ap- 
parent, for example, pertaining to the way in which in- 
dividuals find new job opportunities. It has been consis- 
tently shown by sociologists and economists alike 0, ^ 
that personal acquaintances or neighborhood effects play 
a prominent role in job search. This, in turn, leads to 
significant correlation in employment across friends, rel- 
atives, or neighbors. The common thesis proposed to 
explain this evidence is that, in the presence of environ- 
mental volatility, the quantity and quality of one's social 
links - sometimes referred to as her social capital ,12i]- is 
a key basis for search and adaptability to change. 
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The study of complex networks has been mainly con- 
cerned with simple phenomenological models reproduc- 
ing the main stylized facts. In contrast, the socioeco- 
nomic literature has studied game-theoretic models of 
network formation [l^ in a fixed environment. Here, 
our objective is to integrate and enrich both approaches, 
proposing a stylized model of a society that embodies the 
following three features: (i) agent interaction, (ii) search 
cum adjustment, and (iii) volatility (i.e. random link re- 
moval). Individuals are involved in bilateral interaction, 
as reflected by the prevailing network. Through occa- 
sional update, some of the existing links have their value 
deteriorate and are therefore lost. In contrast, the indi- 
viduals also receive opportunities to search that, when 
successful, allow the establishment of fresh new links. 
Over time, this leads to an evolving social network that 
is always adapting to changing conditions. The model 
aims to capture the continuous struggle of search against 
volatility, which is at the heart of network's dynamics. 
In the long run, the so-called Red Queen Principle Q 
applies: "...it takes all the running you can do, to keep 
in the same place." Agents' continuous search must be 
strong enough to offset volatility if a dense and effective 
social network is to be preserved. On the other hand, 
search can be effective only in a densely networked soci- 
ety. So information diffusion and a dense network of in- 
teractions are two elements of a feedback self-reinforcing 
loop. As a result, the system displays a discontinuous 
phase transition and hysteresis, enjoying some resistance 
to a moderate deterioration of the underlying environ- 
mental conditions. Such a resilience can be interpreted 
as consequence of the buffer effects and enhanced flexibil- 
ity enjoyed by a society that has accumulated high levels 
of social capital. 

These features are captured by a mean field theory 
which is in good agreement with numerical simulation re- 
sults. This theory highlights the particular role that clus- 
tering plays in the dynamics of the model. Indeed search 
is particularly effective when clustering is low whereas it 
is suppressed in a high clustered society. 

The model introduced here is a simplification of a more 



2 



elaborate model proposed by one of us in [T^ to under- 
stand how the network dynamics impinges on strategic 
behavior. The model is also similar to that proposed in 
Ref. 10 to explain the emergence of the small-world 
property in social networks. There are, however, im- 
portant differences between Q and our approach, as we 
shall discuss later at the end of the paper. Related issues, 
within the vast recent literature on network dynamics, 
have also been addressed in Ref. 0|, that studied the 
evolution of network among agents involved in an iter- 
ated Prisoner's Dilemma, and in Ref. that found a 
topological phase transition in networks that minimize a 
suitably chosen cost function. 

Formally, the network is given by a set of nodes TV 
and the corresponding adjacency matrix A with elements 
A,j G {0,1}, for i,j e {1,2, ...,7V}. The value Aij = 1 
means that there is a link connecting nodes i and j , while 
Aij = holds otherwise. We will require An = (no on- 
site loops) and Aij — Aji (unoriented links). The matrix 
stochastic process A{t) in continuous time t represents 
the evolution of the network. 

Two local parameters will be of central importance for 
our discussion namely the node degree Ci (t) — Aij (t) 
and the local clustering coefficient 



q^{t) 
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The latter measures the fraction of pairs of neighbors of i 
who are also neighbors among themselves. The averages 
of these two quantities over sites will be simply denoted c 
and q. While random networks have q 1/iV, social net- 
works typically have a clustering coefhcient Q bounded 
above zero. 

Denote by Fi — {j\Aij — 1} the set of neighbors 
("friends") of the node i. The network evolves due to 
the following three processes. 

1. Long distance search: At rate 77, each node i gets 
the opportunity to make a link to a node j randomly 
selected (if the link is already there nothing happens). 

2. Short distance search: At rate ^, each node i picks at 
random one of its neighbors j € Fi and j then randomly 
selects (i.e. "refers to" ) one of its other neighbors k G 
Fj\{i}. li k ^ Fi then the link between i and k is formed. 
li Fi = or Fj = {i} or k € Fi nothing happens. 

3. Decay: At rate A, each existing link decays and it 
is randomly deleted. 

For ^ = 0, the dynamics is very simple and the sta- 
tionary network is a random graph with average degree 
c = ry/A. For 77 <C A the network is composed of many 
disconnected parts. Fig. ^ reports what happens when 
the local search rate ^ is turned on. For small ^, net- 
work growth is limited by the global search process that 
proceeds at rate rj. Clusters of more than 2 nodes are 
rare and when they form local search quickly saturates 
the possibilities of forming new links. Suddenly, at a 
critical value ^ giant component connecting a finite 
fraction of the nodes emerges. The average degree c in- 
deed jumps abruptly at The distribution p{c) of c; 



is peaked with an exponential decrease for large c and a 
power law p{c) ~ c'^ for c small. The network becomes 
more and more densely connected as ^ increases further. 
But when ^ decreases, we observe that the giant com- 
ponent remains stable also beyond the transition point 
(C < ^2)- Only at a second point £,1 does the network 
lose stability and the population gets back to an uncon- 
nected state. There is a whole interval [■Cij C2] where both 
a dense-network phase and one with a nearly empty net- 
work coexist. This behavior is typical of first-order phase 
transitions. The coexistence region [^1,^2] shrinks as 77 
increases and it disappears for 77 > 0.05A. 
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FIG. 1: Average degree c (top) and clustering coefficient q 
(bottom) from numerical simulations with rj/X = 0.01 for 
populations of size n — 1000 and 2000. Runs were equili- 
brated for a time teq = 3000/ A before taking averages. 



The average clustering coefficient q shows a non-trivial 
behavior. In the unconnected phase, q increases with ^ 
as expected. In this phase, q is close to one because the 
expansion of the network is mostly carried out through 
global search, and local search quickly saturates all pos- 
sibilities of new connections. On the other hand, in the 
dense-network phase, q takes relatively small values. This 
makes local search very effective. Remarkably we find 
that q decreases with ^ in this phase, which is rather 
counterintuitive: increasing the rate £ at which bonds 
between neighbors form through local search, the density 
q of these bonds decreases. In fact, similar behavior is 
found if, fixing ^ and r/, the volatility rate A decreases. 

To shed light on these numerical results, we study the 
dynamics of the distribution P(c, q, t) of the degrees c 
and clustering coefficients q. Specifically, we study a 
mean field approximation that assumes qi — q for all 
i — 1 , . . . , TV and 
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P{c,q,t) = Y[p{c^,t)S{q^ - q). 
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It is convenient to set A = 1, by an appropriate time 
rescaling. Then, the transition rates that enter into the 
master equation for p have the form: 



w(c ■ 



c+ 1) 
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c 



- (3e{c) + 7c 



(3) 
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where e{k) = for fc < and e{k) = 1 for fc > 0. In 
Eq. (PI, the term 2ri accounts for fong distance search. 
The factor 2 counts both the processes when the search 
opportunity is given to site i and when it is given to 
another site, which selects i as end point. The second 
term in Eq. arises from focal search and it requires 
that Ci > 1. Here /3 = ^(1 — q)P{cj > l\j € Fi} accounts 
for the fact that the selected friend j £ Fi must have at 
least one more friend k ^ i and that k ^ Fi, which occurs 
with probability I — q. 

Finally the last term accounts for indirect local search 
opportunities given to a friend fc of a friend j € Fi of i. 
This process is proportional to Ci and 7 = C(l ~ 'l){'^~k^) 
accounts for the probability that i ^ Fk and k selects 
j. Note that the probability that j G F^ selects i is 
l/{cj — 1). This, combined with the multiplicity Cj — 1 
of second neighbors k 6 Fj of i, contributes a factor 
{cj — l)/(cj — 1) = 1 to the rate. 

Both P and 7 will be determined self-consistently. The 
master equation for p{c) , in the stationary state, can be 
solved using the generating function tt{s) — {s'^): 



7r(s) 



/3 + 2t/(1 --75)^^ _ 2r/ + f3 

/? + 277(1 - 7)-f ' ^~ 7 ■ 



(5) 



Note that p{c) ~ for small c and p{c) 



,-\lnj\c fQj. 

large c which perfectly matches the behavior observed in 
numerical simulations. 

Eq. (O, allows us to compute the distribution P{cj = 
k\j £ Fi} = p{k) for the degree Cj of j £ Fi. The larger 
Cj the more likely is j a neighbor of i. Thus, p{k) oc 
kp{k) which, in terms of the generating function, implies 
7r(s) = S7r'(s)/7r'(l). Using this to compute averages over 
Cj and Cfe we arrive at the self-consistent equations: 



q) 



1 



^'(1) 

1 - 7r(0) 



(6) 
(7) 



As we should, in the limit ^ — + we find /3, 7 — > and we 
recover a pure Poisson distribution with mean 277. But 
with constant q Eqs. H()I7|I arc not able to reproduce the 
observed behavior. It just predicts a smooth crossover 
and no phase coexistence. This means that, in order to 
shed light on our observations, it is essential to allow 
for q to depend on the parameters of the model and the 
distribution p{k). 



In order to derive an equation for g, we consider the 
evolution of clustering for a node i with Ci friends. Let 
Qi = J2j<k ^ij^ikAjk = qiCi{ci - l)/2 be the number of 
pairs of friends of i which are also friends among them- 
selves. Only local search processes contributes to an in- 
crease in Qi through two different routes. 

The first is when a local search opportunity is given 
to site i itself and has already been discussed above. 
Its rate is Wi{Q^ ^ + 1) = ^(1 - q)P{c^ > 1}. 
The second occurs when a local search opportunity is 
given to some j £ Fi, who then asks to i about his 
other friends k £ Fi [k ^ j). This may lead to the 
formation of the link between j and k, thus increas- 
ing Qi by one. The rate of this process is given by 
W2{Qr ^ + 1) = i{c,9{c, - l))(c-i)(l - q). Here, 
1/cj is the probability that j picks i from his neighbors 
and 1 — g is the probability that k ^ Fj. This rate 
should is multiplied by the number Ci of neighbors of 
i, but is zero unless Ci > 2 [Ts] . Finally, we must account 
for the link-decay process that, contrary to the former 
two, decreases Qi. The rate at which this happens is 
Wx{Q, ^ - 1) = = (c,(q - 1)) q/2. 

If we now take the averages on Ci and Cj with probabil- 
ity distributions p{k) and p{k) respectively, we can im- 
pose stationarity on Qi, i.e. (AQi) = W1+W2 — Wx ~ 0. 
After some algebra, this condition becomes the desired 
equation for q: 
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e(l-g)[2-40)] 



1 



V(0) 



^'(1) 



(8) 



Using Eq. (|3Jl, we arrive at a set of three self-consistent 
equations for (3, 7 and q. 




FIG. 2: Mean degree c (top) and clustering coefficient q (mid- 
dle) as a function of ^/A for 77/A = 0.19. Bottom: phase di- 
agram in the mean field approximation. The behaviour of c 
and q along the dashed line is reported above. 

The solution of Eqs. ^ 0) and lO is shown in Fig. 
121 for 77 = 0.19. The average degree c = 7r'(l) and clus- 
tering coefficient q exhibit the same dependence on ^ as 
that found in numerical simulations. Depending on the 
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values of ^ and rj we find either a unique solution witii 
small c and large q (corresponding to a dilute network) 
or a unique solution with large c and small q (dense net- 
work), or both solutions simultaneously. In particular 
the solution correctly reproduces the behavior of q in 
the two phases: q increases with ^ in the dilute network 
phase whereas it decreases with ^ when a giant compo- 
nent forms. Our approach shows that this is not just a 
by-product of our analysis but rather an essential ingre- 
dient for understanding the network's dynamics. Fig. |21 
also depicts the phase diagram predicted by the mean 
field. In the shaded region, Eqs. lIHl have three 
solutions, of which two are stable and correspond to the 
two coexisting phases. 

The coexistence interval [■^i, 1^2] predicted by the mean 
field theory is much larger than that observed in numer- 
ical simulations, as can be seen by comparing Fig. ^and 
Fig. 121 For example, the critical point rj^/X = 0.226. . . 
above which there is a smooth crossover in the mean field 
is an order of magnitude larger than that suggested by 
numerical simulations. We believe this is due to the fact 
that mean field theory underestimates fluctuations and 
neglects correlations. 

It is straightforward to repeat our approach to ob- 
tain the mean field equations for the model of Ref. . 
The network growth process of Ref. [Tl mixes local and 
global search in a different way. (For example, search is 
always effective since when it is unsuccessful locally, the 
agent nevertheless creates a link through global search.) 
In addition volatility affects sites instead of removing 
bonds, i.e. Eq. I^J is replaced with w(ci ^ 1) = p. This 
changes considerably the stationary state distribution, 
since for c > 1 the stationarity in the master equation 



implies p{c) = a+p/7+c P('"~' -'^)' ^-^^ ^ POwer law behavior 
p{c) ~ c~^'/'*'~^, as observed in Ref. The solution of 
the self-consistent equations is always unique, implying 
that there no phase transition in this model. These con- 
clusions illustrate that the present approach is a rather 
powerful tool in the analysis of network dynamics in a 
wide range of different setups. 

To sum up, our aim in this paper has been to study 
a simple model of network formation whose implications 
shed some light on the evolution (rise and fall) of a net- 
worked society. The induced network dynamics displays 
a first-order transition that, as the environmental condi- 
tions improve, lead from a sparse phase to a qualitatively 
different regime where the rich potential of a network so- 
ciety is realized. Thus, in particular, social interaction 
becomes dense (i.e. average connectivity is high) and in- 
dividual search turns effective (i.e. redundant search is 
avoided by low clustering). 

These findings explains the apparently paradoxical ob- 
servation that a networked society does not necessar- 
ily materialize even under favorable conditions while, by 
contrast, it displays a significant resilience to deteriorat- 
ing conditions. This may help understand the origin of 
the "miracles" and "anti-miracles" in economic develop- 
ment |20j , which are still an unresolved puzzle for modern 
economic theory. 
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